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Abstract 

We define some approximation schemes for different kinds of generalized backward stochastic 
differential systems, considered in the Markovian framework. We propose a mixed approxima- 
tion scheme for the following backward stochastic variational inequality 

dYt + Fit, Xt,Yt,Zt)dt + Git, Xt,Yt)dAt e d^iYt)dt + ZtdWt, 

where (Xt, At)t^[Q T] is the unique solution of a reflected forward stochastic differential equation. 
More precisely, we use an Euler scheme type for the system of decoupled forward-backward 
variational inequality, combined with Yosida approximation techniques. 

Key words and ptirases: Euler scheme, Yosida approximation, error estimate, multivalued back- 
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1 Introduction 

The stochastic differential equations (SDE) with reflecting boundary conditions, also called reflected 
stoctiastic differential equations (RSDE) appears from the modeling of different kinds of constrained 
phenomenon. The elliptic and parabolic PDEs with Neumann type and mixed boundary conditions 
lead us to probabilistic interpretations, via the Feynman-Kag formula, in terms of reflected diffusion 
processes, which are solutions of RSDEs. This type of equations were studied for the first time by 
Skorohod in [22] and after this, considered in general domains (see [11], [TJ], [21], [23]...). 

Since 1990s a lot of researchers focused their attention to numerical schemes, methods and 
algorithms for the study of the behavior of the solution for RSDEs. In the recent years, some 
new techniques consist in splitting-step algorithms and mixed penalization methods. The Euler 
approximation was considered for the first time by Chitashvili and Lazrieva in [5], followed by 
the Euler- Peano approximation, which was introduced by Saisho in [21j. Lepingle in [10] and 
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Slominski in ^23j analyzed the corresponding numerical schemes and their rates of convergence. In 
order to approximate the solution of RSDEs, the penalization method was also very useful (see 
Menaldi, [H])- Approximation methods for a diffusion reflected and stopped at the boundary 
appear in the literature in 1998, in the paper of Constantini, Pacchiarotti and Sartoretto [6]. They 
defined a standard Euler projected approach to stopped reflected diffusions, approach which yields 
a method with weak order of convergence (1/2 in particular) and they give an easy example where 
this convergence rate is precise. Regarding adaptive approximations of one-dimensional reflected 
Brownian motion, it can be used a simple method of two fixed step sizes chosen according to the 
distance at the boundary. 

In the paper p], Asiminoaei and Ra§canu used a mixed method consisting in penalization and 
splitting-up for the study of multivalued SDE with reflection at the boundary of the domain. The 
penalization method was also used by Ra§canu in [19] for the study of the generalized Skorohod 
problem and of its link to multivalued SDE governed by a general maximal monotone operator (of 
sudifferential type). Recently, in [7], Ding and Zhang combined the penalization technique with 
the splitting-step idea to propose some new schemes for the RSDE in the upper half space. 

In 1990, Pardoux and Peng introduced in [15] the notion of nonlinear backward stochastic 
differential equation (for short, BSDE), and they obtained the existence and uniqueness result for 
this kind of equation. Since then, the interest in BSDEs has kept growing and there have been a lot 
of works on that subject, both in the direction of the generalization of the equations that appear and 
in constructing schemes of approximation for them. The backward stochastic variational inequalities 
(for short, BSVI) were analyzed by Pardoux and Ra§canu in [17J and [18j (the extension for Hilbert 
spaces case) by a method that consists in a penalizing scheme, followed by its convergence. 

Starting with the paper of Pardoux and Peng |16j . have been given a stochastic approach for 
the existence problem of a solution for many type of deterministic partial differential equations 
(PDE for short). In |17j it is proved, using a probabilistic interpretation, the existence for the 
viscosity solution for a multivalued PDE (with sub differential operator) of parabolic and elliptic 
type. More recently, Maticiuc and Ra§canu in [12], prove an extended result concerning generalized 
type of BSDE (including an integral with respect to an adapted continuous increasing function and 
two subdifferential operators). These type of BSVI allows to prove Feynman-Kac type formula 
for the representation of the solution of PVI with mixed nonlinear multivalued Neumann-Dirichlet 
boundary conditions. 

Even this type of the penalization approach is very useful when we deal with multivalued 
backward stochastic dynamical systems governed by a subdifferential operator, it fails for the case 
of a general maximal monotone operator. This motivated a new approach, via convex analysis, for 
the study of both forward and backward multivalued differential systems. In [20J, Ra§canu and 
Rotenstein identified the solutions of those type of equations with the minimum points of some 
proper, convex, lower semicontinuous functions, defined on well chosen Banach spaces. 

Euler-type approximation schemes for BSDE, and for BSDE with exit time for the forward part 
of the system, were introduced by Bouchard and Touzi in [4J and Bouchard and Menozzi in [3]. 
They considered the Markovian framework of a coupled forward-backward stochastic differential 
system and they defined an adapted backward Euler scheme for the strong approximation of the 
backward SDE with finite stopping time horizon, namely the first exit time of the forward SDE from 
a cylindrical domain. In [2j, Bouchard and Chassagneux study the discrete-time approximation of 
the solution of a BSDE with a reflecting barrier. 

The paper is organized as follows. Section 2 presents some basic notations, hypothesis and 



2 



results that are used throughout this paper. Section 3 is dedicated to the analysis of the behavior 
of an approximation scheme defined for a backward stochastic variational inequality. In Section 
4 we present an existence and uniqueness result for a generalized BSVI and we propose a mixed 
Euler type approximation scheme for its solution. 



2 Notations. Hypothesis. Preliminaries 

In all that follows we shall consider a finite horizon T > and a complete probability space 
(fijj-", P) on which is defined a standard d-dimensional Brownian motion W = (Wt)f^rp whose 
natural filtration is denoted F ={J-^t, < t < T}. More precisely, F is the filtration generated by 
the process W and augmented by A/jp, the set of all P-null sets, i.e. Tt = (^{Wg, s <t} V A/p. 

We denote by U^^{Q; C([0, T] ; M'^)), r G [1, oo), the closed linear subspace of adapted stochastic 
processes / G L''{Q,T,F;C{[0,T];R'')), i.e. /(•,*) : 17 ^ is J^-measurable for all t G [0,r] 

and E ^supigjo,r] 1/ (*)r) < Also, we shall use the notation Lj;^(fi; L9([0, T] ; R'^)), r,g G 
[l,oo) the Banach space of J^t-measurable stochastic processes f : n X [0,T] ^ R'' such that 
E(/o^|/(t)rdt)'^''<oo. 

Consider the following data: 

• the continuous coefficient functions 6 : [0, T] x M™ ^ R™, a : [0, T] x ^ ]R™>''^, 5 : M"^ ^ 
R" and F : [0, T] x R™ x R" x R"^'^ — )- R, which satisfies the following standard assumptions: 

for some constants q G R, L, /3, 7 > and for all t £ [0,T] , x, x £ R™', y, y G R" and z, 



z G R"""^ 



(i) \b (t, x) — b (t, x)\ + \\a (t, x) — a {t, x)\\ < L\x — x\ 



(1) (a) {y-y,F{t,x,y,z) -F{t,x,y,z)) <a\y-yf , 

{in) \F{t, X, y, z) - F{t, x,y,z)\ < /3 ||z - z\\ , 

and there exist some constants M > and p, q £ N such that, for all t G [0,T], x G R™" and 
y G R" : 

(2) \g{x)\ < M{l + \x\'^) and \F{t,x,y,0)\ < M{1 + \x\p + \y\). 

• the function ip : R" — )• (—00, +00] which is a proper convex lower semicontinuous function 
and satisfies that there exist M > and r G N such that, for all x G R™" : 

(3) M9{^))\<M{l + \xn. 

The following theorem summarizes some already well known results concerning forward and 
backward SDE, considered in the Markovian framework (for the proof see Karatzas & Shreve [8], 
for forward case, and Pardoux & Ra§canu [17] for the backward system). 

Theorem 1 Let (t, x) G [0, T] x R™ be fixed. Under the assumptions (CP, ^ and the forward- 
backward coupled system 

' dXt'"" = bis,Xt''')ds + , s G [0,r] , 

(4) I dn*'" + F{s, X*'", y,*'", G dip(Ys'''-')ds + Zl^'^dWs , t G [0, T] , 
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has a unique solution, i.e. there exist a unique process X*'^' G L^^{O,;C{[0,T] such that 

(5) X*'^ = a;+/ h[r,Xl^^)dr+ a(r, , s G [0, T] , 

Jt Jt 

and respectively 

(l^*'^^*'^f/*•") e^L(^^;C([0,r];M")) x LUn;L\[0,T];R-'<'')) x L^n; L\[0,T] ■,R^)), 
such that 

yi'^+ / U','^dr = g{X'/)+ f ^,T]{r)F{r,Xl'\Y^^^,Zlndr- j Z^'W, , s e [0, T] 
[7^'"= G 5^9 fy,*'"^) , dPxds on X [0, T] . 



Moreover, for all p > 2, there exists some constant Cp > 0, (7 G N*, such that, for all t, t £ [0,T] , 
X, X e M" : 

0') E(sup,g[o,T]|^s'1^) <Cp(l + |xn, 

O'i) E(sup,g[o,T] l^s'" - ^I'T) < + 1^1'' + - + 1^ - ^1'')' 

(jjj) E(sup,g[o,T] < Cil + \x\'), 

(jv) E( sup,g[o,T] - y/'^n < C2 [^g{X'/) - g{X'/)\^ + 

+E^^|l[,,^](r)F(r,X*•^y,*■^Z*'^) - l[,-^^](r)F(r,Xp,y/'^ V . 

3 Approximations schemes for BSVI 

We will consider a partition of [0, T], 

TT = {ti = ih : <i <n}, with /i := T/n, n G N*, 

on which we approximate the solution of the backward stochastic variational inequality ([6]) . For the 
numerical simulations of the forward part, the most standard approach consists in approximating 
the SDE in a proper way on each interval [ti,ti^i] by the classical Euler scheme (see, e.g. Kloeden 
& Platen [9]): 

f =Xu+b {X,^) h + a {X^) {Wu^, - Wu) , ^ = 

[ X^' = Xq. 

We remark that the above numerical scheme is easy to implement since it requires only the simu- 
lation of d-independent Gaussian variables for the Brownian increments, providing a weak error of 
h order. 

For t G [tj, tj+i] let 

X^ = Xl + b{X^;) {t - t,) + a{Xl) [Wt - Wu) . 
We have the following estimation of the error given by the Euler scheme (see [9]). 
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Proposition 2 Under the assumptions (CP on the coefficients b and a, for all p > 1, there exists 



Cp > such that 



0,n-l Vfe[0,T] 



max El sup \Xt - + sup \Xt - X, 



i+ij 



i/p 



< Cr, Vh. 



Here and subsequently we will consider the one-dimensional BSDE case. 

Using the Yosida approximation Vipe of the multivalued operator dip, with e = h"" and a G 
(0,1/2) (the way of choosing this constant will be detailed later), we deduce that the following 
approximate equation 



(7) ¥1"+ I V^h-{Yr)dr = g{XT)+ [ F {r , Xr , Z^)dr - [ Z^dWr, Vt £ [0,T] 
Jt Jt Jt 



a.s.. 



admits a unique solution {Y/',Zj^) G Ll^{n-C{[0,T];R)) x Ll^{n; L'^{[0,T] ;] 

Further, inspired by the paper of Bouchard and Touzi, ^4], let us define an Euler type approx- 
imation for the Yosida approximation process y^. For an intuitive introduction, let Yp := g{X^) 
be the initial condition, and, for i = n — 1,0, remark that 



(8) 



~ Yt , + h 



F{t,,XlYll,Z^^) - Viph'^{Yll)\ - Z^iWu^, - Wu) 
taking the conditional expectation E* (•) := E (• \J-ti), we obtain 

Yll ~ EX^iJ + h [F{t,,XlYll,Z^^) - Viph^Y/^) 
If we multiply ([8]) by Wt^^-^ — Wt^ it follows 

z^^^lriY,l^{Wu^,-Wu)). 

Therefore, we propose the following implicit discretization procedure, which define the pair (Y^, Z^^ 
inductively, for i = n — 1,0 : 

( Y^ :=g{X^), Z^ = 0, 



(9) 



ti \ ti+i ' 

1, 



F(t„x/^,y,^,z/^)-V(p,.(y,^ 



where E 



i,h I 



Zl\:= -W^\Yll^^{Wu^,-Wt,)), 
y U,]:=Vipf,a{W^HYt,)), 
E(. IJ-/;) and J-/^:=a(X/;, :0<i<i). 



Remark 3 Observe that Y^^ is defined implicitly as the solution of a fixed point problem. Since 
the involved functions are Lipschitz, it is well defined. Moreover, for small values of h > it can 
be estimated numerically in an accurate way. 



5 



Remark 4 We can also use an explicit scheme to define 

The advantage of this scheme is that it does not require a fixed point procedure but, from a numerical 
point of view, adding a term in the conditional expectation makes it more difficult to estimate. 
Therefore the implicit scheme can be more tractable in practice. 

Remark 5 We have that the filtration J^t generated by the Brownian motion coincides with the 
filtration generated by the diffusion process X, i.e. Tt = ^ and, from the Markov property of the 
process X^, it follows that 

W{Y,l^) = W'^{Y,lJ = E{Y,l^ \X^J, 

^\Y,l^iWu+i - Wu)) = E^'HYt,{Wu+, - Wu)) = E{Y,l^iWu^, - W^) \X^). 

Consider now a continuous version of ([9]). From the martingale representation theorem there 
exists a square integrable process such that 



(10) 



Y,l^=E\Y,'^ )+ / Z^dWs, 



ti + l 



and, therefore, we define, for t £ 



(11) 



Yl^:=Y,'l-{t-t, 



f{t,,XlYll,Zi\)-V^h^{Yll) 



+ / zUWs. 



Obviously, we obtain that 1^^ = 1^^, and, for the simplicity of the notation, we will continue to 



write YI' for YI'. 



Remark 6 From MU\) and the isometry property, we notice that, for i = 0,n — 1, 



(12) h zi = r(yi^(vFi,^, - Wt,)) = r 



(m,+i-mj/ zidw, 



To approximate Z^ we use 



z^ := -e:' 



Z^ds 



rather than Z^_, which is the best approximation in ($7 x [0, T]) of Z'^ by adapted processes 
which are constant on each interval [tj,tj+i) (see Lemma 3.4.2 from Zhang [2^): 



E 



ti+i 



\Z^ - Zj^.l'^ds 



< E 


[/ 







'■ti+l 



for all J^i- -measurable stochastic process ry. 
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Remark 7 From (T^) . the definition of Z/* and Jensen inequality we obtain 



\ Li Li \ 



(13) 



ti+l 



E'lA'^Zsl'^ds = f 
h 



E 



-I 2 



A'^Zsds 



< -E 
h 



ti+l 



ds < 



ElA^ZJ'^ds. 



In order to prove an error estimate of the scheme first we use the solution (Y/^, Zj^^^^^^^ of the 
approximating equation ([7]). The next result is a straightforward consequence of Proposition 2.3 
from Pardoux & Rascanu [17j . 

Proposition 8 Under the assumptions ([2P-(0j; there exists C > such that 

E\Yt-Yl'\'^ + E I 

te[o,r] 

rT 



(14) sup EiFt -y^^^l^ +E / \Zt- Z'l\^dt<CT{T)h"', 

;e[o,r] Jo 

where V{T) := E 1 + |c/(Xt)|^ + \Xt\^ + F ds 



We recall Theorem 3.4.3 from [23], applied for the solution {Y^^ , Z^^ of d?]). To otain a similar 
conclusion we have to impose more restrictive assumptions than (JTJlSj): 

• there exists some constant K > such that 

(i) \h{x) -h{x)\ + \\a{x) - a {x)\\ <K\x-x\, Vx,x G M", 



(15) 



{ii) |F(o-F(e)i <K\i-i\, ve,eG [o,t] x 

(m) \g{y)- g{y)\<K\y-y\,\ly,y 



X K X 



• the function ip : — )• (— oo, +oo\ is a proper convex lower semicontinuous function and there 
exist M > and r G N such that 



(16) 



W{g{x))\ < M(l + IxD, Vx G 



Proposition 9 Let the assumptions ( fi5j) and i fifi]) 6e satisfied. We have the following estimate, 
for some C > 0, 



max sup 

=0,n-lt(^[ti,U+ 



where Z^^ := ^E' Z^ds 
Proof. The inequality 



" fU+i 

ElY/" - Y/l^^\^ + ^E / \Z^- Zj^fds < Ch 
i] j=i Jti 



max sup EIF^'^ - F^^^ < 



can be obtained by classical calculus, using Ito's formula, Lipschitz property of the coefficient 
functions and the bounds of the approximate solution (Y'j'^, Z/')^^^ of ([7]) (see Proposition 2.1 
and 2.2 from [I?]). 
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For the proof of the inequality 



- Z^]^ds < Ch 



i=l 

is sufficient to recall the proof of Theorem 3.4.3 from 

Using the estimates from the above Propositions we can prove the following: 

Proposition 10 Let the assumptions ( liJ)) and ( fjfij) be satisfied. Then there exists C > such that 



sup Ely^'* - yI"]^ + E I |Zf - Zl^l'^dt < Ch 

tG[0,T] Jo 

Proof. From ([7]) and (fTTj) we deduce that, for i = 0, n — 1 and t G 



1-2a 



Y/" = Y." + 



ti + 1 



ti+1 



Fis,X,,Y,\Z^)-VipH^{Y: 



F{t„X^.Y/^,Z^)-V^^aiY/^^ 



ti+l 



ds - 
ds 



Z^dWs, 



ti+1 



Z^dWs. 



Throughout the proof let A^Ft := F{t, Xt,Yt'' , Zj^) - F{ti,X^^,Yll, Zl^^), A'^Yt := Yt^ - Y{^ and 
l^Zi := Z^ — Z^, t € 



t 

Applying Energy equality we obtain that 

''ti + l 



(17) ElA^Ftp + I ¥.\A^Zs\^ds = ElA'^Ft^+J^ + 2E / A!'Ys A^F^ds 

Jt Jt 

- 2E^*'''' A'^Ys (Vifh^Yi^) - Viph4yl!)) ds, 

We first compute A'^Ys ■ A^F, - {Viph^Y,^) - ViphaiY/")) for which we use Lipschitz property 

of F and 'Vcph^ : 

(18) 

Ia'^F, - (V^H^Yh - V^h^{Yll)) 



2E 
< 2i^E 



A^y,- 



A^Y. 



ds < 



1 



Is - 1,\ + \Xs -x^^\ + \y:^ - + \z^ - z/^i + ^\Yl^ - yH^ 

'■ti+1 



ds < 



< (K^a + I3)E [ \A^'Ys\^ds + -E [ - Pds + -h^+ 

Jt Jt 'a 



-1e 

a 



ti+1 /A 1 

\Z^-Zt\^ds+(- + -^ 



E 



ti+1 



\Y,'-Yl^\'ds, 



where a, f3 > will be chosen later. 

From now on, let C > be a constant independent of h, constant which can take different 
values from one line to another. 

From Proposition [2] we have that there exists C > such that 



E\Xs - Xj^]'^ < 2E\Xs - XtS" + 2E\Xt^ - X]^]' < Ch, 



h\2 
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and, from Proposition [U 

E\Y,^ - yII\^ < 2^Y^ - Vl'f + 2^yII - yII\^ <Ch + 2^/^^Yt^\^. 

Using ^ 



¥.\Z^ - Z^f < 2E\Z^ - Z^f + 2E|Zj'^ - Zji'f = 2E\Z^ - Z^f j EjA'^Z^pds. 

J to 



Then ([T8]) yields 

(19) Ai {t) := E\A^Yt\^ + / E\A''Zs\^ds < {K^a + /3) E / \A^Y,\^ds + Bi 

Jt Jt 



where 



(20) 



:= E\A''Yt^_^f + - E\Z^' - Z^l^ds + - T^' EjA'^Z.^ds + f- + 



+2/i(l + ^)E|A%l' 



Using a backward Gronwah type inequality we deduce 



ElA'^y^P < Bie 



{K^a + I3)ds 



and, therefore, 

(t) <Bi + {K^a + P) CB.ds = B,[l + C {K^a + h] < B, [1 + Ch], he (0, 1) 



The above inequality and the definition of Bi implies 



< (1 + Ch) 



a 



E|A%|2 + 

ti+l 



1 - (1 + C/i) - 



a 



E\A^Zs\'^ds < 



E|A% 1^ + - / E\Z^ - ZlYds + Ch' 



U 



Taking a G (0, 1/2), we can chose h > sufficiently small and a,/3 > large enough such that 



Ci := 1 - (1 + C/i) ( ^ + ) 2h > and := 1 - (1 + C/i) f > 



and, therefore, 



(21) 



CiE\A^Ytf + C2 [ ElA^Z.pds < 
Jt, 

< {1 + Ch) 



E|A"yt,_^J^ + - / E\Z'^ - Z[\\'ds + C/i 



Writing the above inequality for each i = 0,n — 1, we can deduce 

n „+. 
/ ^1- 



" fti+i 



, i = 0, n — 1. 



From the Lipschitz property of g and Proposition [9] we obtain, for each i = 0,n — 1, since a G 
(0,1/2), that 



(22) 



E\A''Ytf < Ch, V/i G (0,1) smaU enough. 



For the proof of the inequahty concerning || A^'Zs ||j^2(Qx[o T]) following manner 

(see, e.g. Bouchard & Touzi [Ij). From ([?T|) it follows that 



ra-l 



(is 



Jo Jt. 

< (1 + Ch) ^E|A'^yj^^j2 ^ c7/j2-2a r EIZ,'^ - Z/; 

= (1 + C7/i) ElA'^yrP + C/i^^^a ^ « / ' g,^/. _ 2^ 

i=o '' 

n-l 

+ ((1 + Ch) + (1 + C/i) (I + ^)2h - l) ^ElA^^y^J^ _ CiElA'^yol 



, n-l 



n-l 



Ci^ElA^y^, 



j=0 



+ 



and, therefore, from ()22p . 

JO 



n-l 



a p a p 



n-l 

^ElA'^YtJ^ < 



-2a 



+ 



1=1 



< c 



n—1 fti-^i 

E\A''Yt\^ + Y. I - Ztfds + /i^-^a ^ ^1- 



2a 



C/i n 



l-2a 



= C[Ch + Ch + Ch^'^"] < Ch 

Using the definition (j20p of we deduce that Bi < Ch and, respectively, maxElA'^Y't-p < Ch, 
which completes the proof. ■ 

Consequently we have proved our main result: 

Theorem 11 There exists the constant C > which depends only on the Lipschitz constants of 
the coefficients, such that: 



(23) 



sup E\Yt - y^^p + E / \Yt- Yl'\' ^\Zt-Zi 
iG[o,r] 



rh\2 



h\2 



dt < Ch^^^^-^^l 
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4 Generalized BSVI. A proposed scheme for numerical approxi- 
mation 

Let 2? be a open bounded subset of M/^ of the form 

V = {x£R'^:£{x)< 0}, Bd {V) = {x G M'^ : £ (x) = 0}, 

where i £ {^'^), |V£(x)| = 1, for ah x G Bd{V). We know that (see, e.g., Lions & Sznitman 
[TT]). for every {t,x) G M+xD, there exists a unique pair of PxM_|_— valued progressively measurable 
continuous processes {Xl'^ , A^s^)s>o, solution of the reflected SDE: 



r-syt psVt i-sVt 

Jt Jt Jt 



(24) 



s I — > AI'^ is increasing. 



A 



t,x 



t,x 



Moreover, it can be proved that 



E sup f - Xf '^'1^) < C{\x - x'f +\t- t'' 
se[o,T] ^ ^ 



At X 

and, for all fi > 0, E(e'^^T ) < oo. 

Consider now the following generalized backward stochastic variational inequality: 

r dYt + F (t, Xt, Yt, Zt) dt + G (t, Xt, Yt) dAt G {Yt) dt + ZtdWt, 0<t<T, 

\YT = g (Xt) • 



(25) 



We will suppose that the functions F and G satisfy the same assumption as the generator function 
F from the previous section. It is known (see Maticiuc & Ra§canu [13j) that the above equation 
admits a unique solution, i.e., for all t G [0,T], P-a.s., 

Yt+ I Usds = g{XT)+ f F{s,X„Ys,Zs)ds+ [ G {s, Xs,Ys) dA^ - [ Z^dWs, 
Jt Jt Jt Jt 

where Ut G dip (Yt) , a.e. on 1^ x [0, T] . 

Theorem 12 Under the considered assumptions, the generalized BSVI 125\) admits a unique so- 
lution (Yt,Zt,Ut) of Ffprogressively measurable processes. Moreover, for any 0<s<t<T,we 
have, for some positive constant C: 



(a) E 



(26) 



where 



j {\Yr\^ + \\Zr\\^)dr + j \Yr\^dAr 
J s J s 



(b) E(^{Yt)) < CM2 and E 



\Ur\ dr 



+ E sup lYrf < GMi, 

s<r<t 



< CM2, 



Ml =E 



leiV J (|F(s,0,0)fds + |G(s,0)fdA 



an. 



d M2 = E(|C|2 + (^(0)- 
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For the generalized system considered above, we propose a mixed approximation scheme, 
considering, for the simpHcity of the presentation, only the case ip = 0. Consider the grid of 
[0,T] : IT = {ti = ih, i < n}, with h := T/n, n G N* and we will define X'^ , the approximating 
Euler scheme for the reflected process X. We follow the paper |6], where the authors present the 
standard projected Euler approach to stopped reflected diffusions. 



xs = 



X, Aq — 0, 

= XI + b{U,X[J{U+, - u) + a{u,xrjiwt,^. 

Taking the projection on the domain, we define 
/ Ki ' Ki e ^> 



Wt. 



and 



"^k+i 



^k ' 



To define an approximation scheme for the generalized BSVI 



consider := g{X^) and. 



for i 



ti+1 



n — 1,0, in an intuitive manner, using the notation AA^. := A^.^_^ — AJ. and AWt^ 

We take the conditional expectation E"^' and it follows 

Yt^ ~ E-^> {Yu^, ) - E-^' [G{Xl^^ , Yu_^,)^Al]. 
This suggest to define the following approximation scheme: 



(27) 



Y- := E^»[y-^^ - G(X-^^,y,,,JAA- 
Zl := ^E^k[Y,^^^AWu - GiX^^^^,Yt^^,)AAlAWu] 



Problem 13 The proof of the convergence for the scheme defined by ^27\ ) can provide a useful tool 
for the approximation of the solution for the Generalized BSVI 125\) . For the moment this is still 
an open problem, which the interested reader is kindly invited to approach it. 
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